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Abstract Neutrino-matter interaction rates are central to the core collapse phe¬ 
nomenon and, perhaps, to the viability of the mechanism of core-collapse 
supernova explosions. In this paper we catalog and discuss the major 
neutrino scattering, absorption, and production processes that together 
influence the outcome of core collapse and the cooling of protoneutron 
stars. These are the essential inputs into the codes used to simulate 
the supernova phenomenon and an understanding of these processes is 
a prerequisite to continuing progress in supernova theory. 
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1. Introduction 

One of the key insights of the 20th Century was that most of the 
elements of nature are created by nuclear processes in stars. Supernova 
explosions are one major means by which these elements are injected 
into the interstellar medium and, hence, into subsequent generations of 
stars. Therefore, supernovae are central to the chemical evolution and 
progressive enrichment of the universe. Core-collapse supernova explo- 
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sions signal the death of a massive star and are some of the most majestic 
and awe-inspiring events in the cosmos. However, to fully understand 
the role of supernovae in the grand synthesis of creation, one must have a 
firm handle on the nuclear data. In particular, understanding the origin 
of iron-peak, r-process, and rp-process elements hinges upon improve¬ 
ments in our knowledge of the properties of exotic nuclei, some of which 
are far from the valley of beta stability. 

The mechanism of core-collapse supernovae is thought to depend upon 
the transfer of energy from the core to the mantle of the inner regions 
of a massive star after it becomes unstable to collapse. Neutrinos seem 
to be the mediators of this energy transfer. In order to understand this 
coupling and the role of neutrinos in supernova explosions, one needs 
to master the particulars of the neutrino-matter scattering, production, 
and absorption rates. Since recently there has been some progress in 
understanding the associated microphysics, it is fitting to summarize 
the neutrino-matter cross sections and the production rates of neutrinos 
in the core-collapse context. To this end, we have assembled here a short 
precis of many of the relevant processes and physics. This contribution 
does not attempt to explain the hydrodynamics of supernova explosions, 
but does try to present the relevant neutrino processes that play a role. 
For the former, the reader is referred to Burrows, Hayes, and Fryxell 
(1995) and Burrows (2000). 

In §2, we present a physical derivation of stimulated absorption (the 
Fermionic correlate to stimulated emission) and then in §3 we present the 
basic cross sections. In §4, we discuss the inelastic neutrino-electron and 
neutrino-nucleon scattering processes and energy redistribution. This 
is followed with a discussion of the alternate, more powerful formalism, 
for determining differential interaction rates in the many-body context 
and for handling redistribution, namely that of dynamical structure fac¬ 
tors. Both the non-interacting and the interacting (in the context of 
a simple nuclear model) cases are explored, as well as collective exci¬ 
tations of the medium. Source terms for electron-positron annihilation 
(§6), neutrino-anti-neutrino annihilation (§7; Janka 1991), and nucleon- 
nucleon bremsstrahlung (§8) cap off our review of the major processes of 
relevance in core-collapse simulations. Nucleon-nucleon bremsstrahlung 
can compete with pair annihilation as a source for i>^, and Vj- 
neutrinos. Clearly, a mastery of these neutrino-matter processes is a 
prerequisite for progress in supernova theory and it is in that spirit that 
we have assembled this review. 
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2. Stimulated Absorption 

The concept of stimulated emission for photons is well understood 
and studied, but the corresponding concept of stimulated absorption for 
neutrinos is not so well appreciated. This may be because its simple 
origin in Fermi blocking and the Pauli exclusion principle in the context 
of net emission is not often explained. The net emission of a neutrino is 
simply the difference between the emissivity and the absorption of the 
medium: 

J^net — • ( 1 ) 

All absorption processes involving fermions will be inhibited by Pauli 
blocking due to final-state occupancy. Hence, rj^ in eq. (1) includes a 
blocking term, (1 — .7v)(Bruenn 1985). J-y is the invariant distribution 
function for the neutrino, whether or not it is in chemical equilibrium. 

We can derive stimulated absorption using Fermi’s Golden rule. For 
example, the net collision term for the process, Ven e-j- e“p, is: 



where p is a four-vector and 

^ e~p) = • (3) 


The final-state blocking terms in eq. (3) are manifest, in particular that 
for the t'e neutrino. Algebraic manipulations convert 'r.[ve.n ^ e~p) in 
eq. (3) into: 

e~p) = 




jFn(l - •^e)(l - ^p) 


J- ^ 




(1 - - .F,, 


( 4 ) 


where 





( 5 ) 


is an equilibrium distribution function for the neutrino and it has 
been assumed that only the electron, proton, and neutron are in thermal 
equilibrium. (I is the difference between the neutron and the proton 
chemical potentials. Note that in there is no explicit reference to a 
neutrino chemical potential, though of course in beta equilibrium it is 
equal to pe — R- There is no need to construct or refer to a neutrino 
chemical potential in neutrino transfer. 
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We see that eq. (4) naturally leads to: 

Jnet=Y^^{B,-L) = n:{B,-I,) . ( 6 ) 

Of course, Bi, is the black body function for neutrinos. This expres¬ 
sion emphasizes the fact that Jnet are the same entity. If 

neutrinos were bosons, we would have found a (1 -|- in the denomi¬ 
nator, but the form of eq. (6) in which Ii, is manifestly driven to B^, the 
equilibrium intensity, would have been retained. From eqs. (4) and (6), 
we see that the stimulated absorption correction to Ka is 1/(1 — 

By writing the collision term in the form of eq. (6), with Ka corrected 
for stimulated absorption, we have a net source term that clearly drives 
lu to equilibrium. The timescale is 1/ ck* . Though the derivation of 
the stimulated absorption correction we have provided here is for the 
i^eU, e~p process, this correction is quite general and applies to all 
neutrino absorption opacities. 

Kirchhoff’s Law, expressing detailed balance, is: 

Ka = rj^lBy or K*a = r\ylB^ , (7) 


where is not corrected for final-state neutrino blocking. Furthermore, 
the net emissivity can be written as the sum of its spontaneous and 
induced components: 


Vv = Ka 


B„ 








( 8 ) 


where -|- or — is used for bosons or fermions, respectively. Eq. (7) can be 
used to convert the absorption cross sections described in §3 into source 
terms. 


3. Neutrino Cross Sections 

Neutrino-matter cross sections, both for scattering and for absorp¬ 
tion, play the central role in neutrino transport. The major processes 
are the super-allowed charged-current absorptions of Ve and Ve. neu¬ 
trinos on free nucleons, neutral-current scattering off of free nucleons, 
alpha particles, and nuclei (Freedman 1974), neutrino-electron/positron 
scattering, neutrino-nucleus absorption, neutrino-neutrino scattering, 
neutrino-antineutrino absorption, and the inverses of various neutrino 
production processes such as nucleon-nucleon bremsstrahlung and the 
modihed URCA process {ve + n + n —> e~ + p + n). Compared with 
photon-matter interactions, neutrino- matter interactions are relatively 
simple functions of incident neutrino energy. Resonances play little or 
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no role and continuum processes dominate. Nice summaries of the 
various neutrino cross sections of relevance in supernova theory are 
given in Tubbs & Schramm (1975) and in Bruenn (1985). In partic¬ 
ular, Bruenn (1985) discusses in detail neutrino-electron scattering and 
neutrino-antineutrino processes using the full energy redistribution for¬ 
malism. He also provides a serviceable approximation to the neutrino- 
nucleus absorption cross section (Fuller 1980; Fuller, Fowler, & Newman 
1982; Aufderheide et al. 1994). For a neutrino energy of ~10 MeV the 
ratio of the charged-current cross section to the i/g“Electron scatter¬ 
ing cross section is ~100. However, neutrino-electron scattering does 
play a role, along with neutrino-nucleon scattering and nucleon-nucleon 
bremsstrahlung, in the energy equilibration of emergent neutrinos 
(Thompson, Burrows, & Horvath 2000). 

Below, we list and discuss many of the absorption and elastic scat¬ 
tering cross sections one needs in detailed supernova calculations. In 
§4 and §5, we provide some straightforward formulae that can be used 
to properly handle inelastic scattering. The set of these processes com¬ 
prises the essential microphysical package for the simulation of neutrino 
atmospheres and core-collapse supernovae. 

3.1 Ue n —> e~ p: 

The cross section per baryon for Vg neutrino absorption on free neu¬ 
trons is larger than that for any other process. Given the large abun¬ 
dance of free neutrons in protoneutron star atmospheres, this process is 
central to Ve neutrino transport. A convenient reference neutrino cross 
section is cTo, given by 

= 1.705 xl0-«cm^ (9) 

7r{ncp 

The total Ve — n absorption cross section is then given by 



( 10 ) 

where gA is the axial-vector coupling constant (~ —1.26), Anp = mnC^ — 
mpCp = 1.29332 MeV, and for a collision in which the electron gets 
all of the kinetic energy Cg- = -|- A^p. Wm is the correction for 

weak magnetism and recoil (Vogel 1984) and is approximately equal to 
(1 -|- l-lEu^/ninC^)- At Eu^ = 20 MeV, this correction is only ~ 2.5%. We 
include it here for symmetry’s sake, since the corresponding correction 
(Wj^) for i?e neutrino absorption on protons is (1 — 7.lEp^/mnC^), which 
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at 20 MeV is a large —15%. To calculate k*, must be multiplied 
by the stimulated absorption correction, 1/(1 — and final-state 

blocking by the electrons and the protons a la eq. (4) must be included. 

3.2 l/g + p —> e+ + n: 

The total Ve — P absorption cross section is given by 


( 1 + \ / Anp \ 

\ 4 y V ruec^ ) 



where e^+ = ey^ — Anp and is the weak magnetism/recoil correc¬ 
tion given in §3.1. Note that is as large as many other corrections 
and should not be ignored. To calculate k*, ap^p must also be corrected 
for stimulated absorption and final-state blocking. However, the sign of 
“ A ™ the stimulated absorption correction for i?e neutrinos is flipped, 
as is the sign of fie in the positron blocking term. Hence, as a consequence 
of the severe electron lepton asymmetry in core-collapse supernovae, 
both coefficients are very close to one. Note that the Pe + P ^ e'^ + n 
process dominates the supernova neutrino signal in proton-rich under¬ 
ground neutrino telescopes on Earth, such as Super Kamiokande, LVD, 
and MACRO, a fact that emphasizes the interesting complementarities 
between emission at the supernova and detection in Cerenkov and scin¬ 
tillation facilities. 


3.3 i^eA A'e 

/,From Bruenn (1985) the total Ve — A absorption cross section, is 
approximated by 



( 12 ) 

where Wuock = (1 - /e-) , Q' = Ma' - Ma + A fin - 

fip + A, A is the energy of the neutron I/ 5/2 state above the ground 
state and is taken to be 3 MeV (Fuller 1982), and the quantities Np{Z) 
and Nn{N) are approximated by: Np{Z) = 0, Z — 20, and 8 for Z < 20, 
20 < Z < 28, and Z > 28, respectively. Nn{N) = 6, 40 — N, and 0 for 
V < 34, 43 < V < 40, and N > 40, respectively. The opacity, corrected 
for stimulated absorption, is then 

K: = XHpNAa%{l-J^n‘!)~^. (13) 



Neutrino-Matter Interaction Rates in Supernovae 


7 


Since Nn{N) = 0 for N > 40, this absorption and emission process plays 
a role only during the very early phase of collapse. Typically at densities 
near p ~ 10^^ g cm“^ k* 0. 

3.4 Wi + p ^ v>i + p: 

The total Vi — p scattering cross section for all neutrino species is: 

(^4 sin^ Ow - 2 Sin^ 0w + j , (14) 

where 9w is the Weinberg angle and sin^ 6w — 0.23. In terms of Cy = 
1/2 + 2sin^ 6w and = 1/2, eq. (14) becomes (Schinder 1990); 

Up = ^ (Cy - if + 3g\{C'^ - 1)^ . (15) 

The differential cross section is: 

where 

p (c(.-i)2 + 3ffi(c;,-i)2 • ^ ^ 

Note that 5p, and 5n below, are negative {5p ~ —0.2 and 5n ~ —0.1) 
and, hence, that these processes are backward-peaked. 

The transport (or momentum-transfer) cross section is simply 

af = ^-^(j^\\{C'y-lf + hg\{C'^-lf] . (18) 

^ D \meC‘^ J L J 

where 


3.5 v>i -\- n —> Ui n: 

The total vi — n scattering cross section is: 

( 1 + A ^ 

4 VrUeC^y y 4 y 


The corresponding differential cross section is: 


dan 

dO. 



( 20 ) 


( 21 ) 



where 


Sn 

The transport cross section is 


9a 

1 + ^9 A 


tr l + 

” 4 \meC^ ) y 6 J 


( 22 ) 


(23) 


The fact that Sp and Sn are negative and, as a consequence, that af' is 
greater than cjj increases the neutrino-matter energy coupling rate for a 
given neutrino flux in the semi-transparent region. 

Horowitz (2002) has recently derived expressions that include a weak 
magnetism/recoil correction analogous to those previously discussed for 
the charged-current absorption rates jZgU ^ pe~ and v'eTi ne"*". We 
take the following form for the weak magnetism/recoil correction, a fit 
to the actual correction factor for the transport cross sections; 

'^n,p ^ T CwM^i^/'^n,p)i (24) 

where, for neutrino-neutron scattering Cwm — —0.766, for neutrino- 
proton scattering Cwm — —1-524, for anti-neutrino-neutron scatter¬ 
ing Cwm — —7.3656, and for anti-neutrino-proton scattering Cwm — 
-6.874. 

In fact, neutrino-nucleon scattering is slightly inelastic and when this 
is germane, as with mu and tau neutrinos, the more general formalisms 
of §4.2 and §5 are necessary. 


3.6 i^i + A —> Ui A: 

In the post-bounce phase, nuclei exist in the unshocked region exte¬ 
rior to the shock. At the high entropies in shocked protoneutron star 
atmospheres there are very few nuclei. There are alpha particles, but 
their fractional abundances are generally low, growing to interesting lev¬ 
els due to reassociation of free nucleons just interior to the shock only at 
late times. However, nuclei predominate on infall and neutrino-nucleus 
scattering (Freedman 1974) is the most important process during the 
lepton trapping phase. 

The differential i^i — A neutral-current scattering cross section may 
be expressed as: 

^ {W Cff + ClosV {•Sion) (1 + 1^) , (25) 

dll b47r \meC‘^J 


W = 1-“ 2sin^ 9w ), 


where 


(26) 
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Z is the atomic number, A is the atomic weight, and (Sion) is the ion- 
ion correlation function, determined mostly by the Coulomb interaction 
between the nuclei during infall. (Sion), in eq. (25) was investigated by 
Horowitz (1997) who approximated it with the expansion 


(‘5ion(c)) 


1 + exp 


-1 -1 


i=0 


(27) 


where 


jZer 1 

a kT 




1/3 


(28) 


a is the interparticle spacing, riion is the number density of ions, T is the 
ratio of Coulomb potential between ions to the thermal energy in the 
medium, and f3i are specified functions of T for each neutrino species. 

Leinson et al. (1988) have investigated the electron polarization cor¬ 
rection, CloSi nnd find that 


^ Z /1 -I- 4 sin^ 6w \ 

= I ( 1 + {kro? ) 

where the Debye radius is 


rn 


I 

AappEp 


(29) 


(30) 


/c^ = |p —= 2(ei,/c)^(l — p), pp and Ep are the electron Fermi 
momentum and energy, and a is the fine-structure constant (~ 137“^). 
Note that rp ^ lOh/pp in the ultra-relativistic limit {pp » mec). The 
Cpos term is important only for low neutrino energies, generally below 
~ 5 MeV. 

Following Tubbs & Schramm (1975) and Burrows, Mazurek, & Lat- 
timer (1981), the form factor term, Cpp, in eq. (25) can be approximated 
by: 

Cpp = , (31) 


where 


y = ^el{r^)l{ncf 




56 MeV 


2 / H \ 2/3 

Viooj 


and (r^)^/^ is the rms radius of the nucleus. Cpp differs from 1 for large 
A and when the de Broglie wavelength of the neutrino is smaller than 
the nuclear radius. 
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When (Sion) = Cff = Clos + 1 = 1, we have simple coherent Freed¬ 
man scattering. The physics of the polarization, ion-ion correlation, and 
form factor corrections to coherent scattering is interesting in its own 
right, but has little effect on supernovae (Bruenn & Mezzacappa 1997). 
The total and transport scattering cross sections for Vi — a scattering 
[Z = 2]A = 4) are simply 

4. Inelastic Neutrino Scattering 

Many authors have studied inelastic neutrino-electron scattering as 
an important energy redistribution process which helps to thermalize 
neutrinos and increase their energetic coupling to matter in supernova 
explosions (Bruenn 1985; Mezzacappa & Bruenn 1993abc). Compara¬ 
tively little attention has been paid to inelastic neutrino-nucleon scatter¬ 
ing. Thompson, Burrows, & Horvath (2000) and Raffelt (2001) showed 
that, at least for mu and tau neutrinos, this process cannot be ignored. 
Here, we review the Legendre expansion formalism for approximating 
the angular dependence of the scattering kernel, detail our own imple¬ 
mentation of scattering terms in the Boltzmann equation, and include a 
discussion of neutrino-nucleon energy redistribution. In §5, we present 
an alternate approach involving dynamical structure factors that is more 
easily generalized to include many-body effects. 

The general collision integral for inelastic scattering may be written 
as 

= (1 - fu) J (33) 

I c(27rhc)3 cos e) 

(34) 

_ -scatt -scatt ^ 

— Hi/ Xu 

where cos0 is the cosine of the scattering angle, Sy is the incident neu¬ 
trino energy, and is the scattered neutrino energy. Although we sup¬ 
press it here, the incident and scattered neutrino phase space distribu¬ 
tion functions {fy and /(,, respectively) have the following dependencies; 
fiy = and /(, = j^y)- ^ and are the cosines of 

the angular coordinate of the zenith angle in spherical symmetry and 
are related to cos 9 through 

cos 9 = ^n' + [(1 — cos((/> — (p'). 



(32) 


(36) 
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The only difference between fi, and J-iy in §2 is that here fi, has explicit 
and El, dependencies, i?™ is the scattering kernel for scattering into 
the bin n) from any bin fi') and is the scattering kernel 
for scattering out of the bin fi) to any bin fi'). The kernels are 
Green functions that connect points in energy and momentum space. 
One may also write R{£y,e'^,cos9) as R{q,(jj), where uj{= — e'^) is 
the energy transfer and q{= + e'^^ — 2eue'^, cos 6]^^^) is the momentum 

transfer, so that the kernel explicitly reflects these dependencies (§5). 

An important simplification comes from detailed balance, a conse¬ 
quence of the fact that these scattering rates must drive the distribution 
to equilibrium. One obtains; i?™ = e~^‘^R°'^^, where (3 = 1/T. There¬ 
fore, we need deal only with The scattering kernels for inelastic 

neutrino-nucleon and neutrino-electron scattering depend in a compli¬ 
cated fashion on scattering angle. For this reason, one generally approx¬ 
imates the angular dependence of the scattering kernel with a truncated 
Legendre series (Bruenn 1985). We take 

°° 2/ -I- 1 

i?°''*(ei.,e'^,cos6») = ^cos 6')Pz(cos 0), (37) 

1=0 ^ 

where 


T + 1 

<h;(e,y,e(,) = J d{cos9) R°'^^{e^,e'j^,cos9)Pi{cos9). (38) 

In practice, one expands only to first order so that 

R°'"\eu,e'^,cos9) ~ ^^oi£u,el) -h ^d>i(ei.,e(,)cos6». (39) 

Substituting into the first term on the right-hand-side of eq. (34) (the 
source) gives 


-scan ^ f 

Jo 


oo riN p' 2 


0 c{2'aTic)^ 


r+1 


1-1 


r27T 

dfi'fl / d^' 

Jo 


1 3 

-<1>0 -h -d>lCOS0 


Substituting for cos 9 using eq. (36) and using the definitions 

1 r+^ 


J. = - 


2 7-1 


dfifu 


(40) 


(41) 


and 


H,. = -r 


1 /•+1 


2 7-1 


dfififu 


we have that 

= (!-/. 


-scatt 


dvr 


c(27r^c)3 7o 


dE'E^Je 




(42) 

(43) 
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Integrating over fx to get the source for the zeroth moment of the trans¬ 
port equation, 


2L 




Similarly, we may write the sink term of the Boltzmann equation collision 
term (second term in eq. 34), employing the Legendre expansion 


-scatt _ 

X.I/ 


deWJ -<ho(l - Jl) - 


ci2TThcyJo " " [2 2 

The contribution to the zeroth moment equation is then 


1 /•+1 


dfii-rr^'fu) = - 


c{2ttTic)^ Jo 


dsyj . 


Combining these equations, we find that 

2 ]-I ^ c{2Tincf Jo " " 


X |i 4>0 [J'(l - J.)e-d^ - (1 - - 1 )| 

(47) 

One can see immediately that including another term in the Legendre ex¬ 
pansion (taking ~ (l/2)4>o + (3/2)$i cos6*-|-(5/2)<h2(l/2)(3cos^ 0 — 
1)) necessitates including P^, and Pi, the second angular moment of the 
neutrino phase-space distribution function, in the source and sink terms. 
While easily doable, we advocate retaining only the linear term and ex¬ 
plore this approximation in the next two subsections. 


4.1 Neutrino-Electron Scattering 

The opacity due to neutrino-electron scattering can be large compared 
with that of other processes at low neutrino energies {su ^ 5 MeV) 
and at high matter temperatures. A good approximation to the total 
scattering cross section has been derived by Bowers & Wilson (1982), 
which interpolates between analytic limits derived in Tubbs &: Schramm 
(1975): 

a, = ^ao{meC^)-^ (Cy + Ca? + J{Cv - Ca? , (48) 

where Cy = l/2-|-2sin^ Ow for electron types, Cy = — l/2-|-2sin^ 6w for 
mu and tau neutrino types, Ca = +1/2 for z/g and and Ca = —1/2 
for Pe and 
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However, the use of cross section (48) implicitly ignores the ineleast- 
icity of neutrino-electron scattering. When inelasticity is germane and 
a full energy redistribution formalism is needed, the scattering kernel 
approach of §4 in which different energy groups are coupled must be 
employed. The scattering kernel R{s^,e'j^,cos6) in §4 is related to the 
fully relativistic structure function for neutrino-electron scattering; 

fl2 

R°'^^{£u,el,cose) = 2G‘^-^[ASi{q,u})+S2{q,u})+BS3{q,uj)]{l-e \ 

(49) 

where (= {u!,q)) is the four-momentum transfer, A = (4e;ye(,-|-g^)/2g^, 

B = Eu + and qa = {£u,qu)- The relativistic structure functions (Si) 
are written in terms of the imaginary parts of the retarded polarization 
functions (Reddy, Prakash, & Lattimer 1998; Thompson, Burrows, Sz 
Horvath 2000; §5): 

5i((7,ti;) = (V^ -I- ^^) ImHf (g, w)-|-Imn-p(g,(j) , (50) 

S2{q,oo) = + J^) lmn.^{q,uj) — J?lmR^{q,u)), (51) 

and 

S2.{q,io) = 2VL4 Imn(^^(g,w). (52) 

V and A are the appropriate vector and axial-vector coupling constants. 
Each of the polarization functions can be written in terms of one-dimensional 
integrals over electron energy {£e), which we label In (Reddy, Prakash, 

& Lattimer 1998) 

2 r 2 " 

lmU^{q,co) = I^+cu/i + ^/o , (53) 

Imn^(g,u;) = l2+‘^h + ^ 

2 

Imn^(g,w) = ^^^0, (55) 

and 

fl2 

Imn(^^(g,w) = [cu/o-h 2/i] . (56) 

Reddy et al. (1998) were able to express the /„’s in terms of polyloga- 
rithmic integrals such that 

/„ = r.(i-h) 


( 57 ) 
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and 


h 



6 



h 





2 ^- 2 ^ + 

z Tz 



(58) 


(59) 


where rje = iXe/T is the electron degeneracy, z = jSuj, uo is the energy 
transfer, and 


In eqs. 
grals; Cn 


OJ 




e_ =-+ - Jl -4^ 

2 2V 


(60) 


(57-59), the ^nS are differences between polylogarithmic inte- 
= Li„(-ai) - Li„(-a 2 ), where 


Un{y) = j dx, (61) 

JO X 

and Lii(x) = ln(l — x). The arguments necessary for computing the 
integrals are ai = exp[/3(e_ -|- cj) — rje] and a 2 = exp(/3e_ — rje). Tables 
for computation of Li„(y) and the I^s can be provided by Sanjay Reddy. 

Figure (1) shows the full scattering kernel for Sy = 20 MeV and = 2, 
10, and 16 MeV as a function of the cosine of the scattering angle, cosO. 
Note that although the absolute value of the energy transfer (|ei, — e(,|) is 
the same for both = 16 MeV and e'^ = 24 MeV, the absolute value of 
i?°"*(20,16, cos 6) is more than twice that of ii°'^*(20, 24, cos 0), reflecting 
the fact that at this temperature the incoming neutrino is more likely 
to downscatter than upscatter. Figure (2) shows the scattering kernel 
for the same conditions as Fig. (1), but also includes both the first-order 
(short dashed lines) and second-order (long dashed lines) approximations 
to We generally employ the former. The latter is included to 

illustrate the improvement in including higher-order terms. In fact, the 
actual degree of expansion necessary to capture accurately the physics 
can only be ascertained by running full transport calculations. We have 
run dynamical simulations with only the zeroth-order and first-order 
terms in the Legendre expansion and find little or no difference between 
the emergent spectra and detailed thermodynamical evolution in the 
models we have studied. Smit (1998) and Smit & Cernohorsky (1995) 
have explored the importance of including the second-order term (oc 
cos^ 6, shown here) and find it negligible. The scattering-angle-averaged 
kernel, also the zeroth-order term in the Legendre series for is 

shown in Fig. (3) for z/g—electron scattering for a matter temperature 
(r) of 6 MeV and with an electron degeneracy factor rje = ye/T = 20 
as a function of for various incoming neutrino energies, e^s. 
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In a full time-dependent simulation it is numerically costly to com¬ 
pute the Legendre moments of the scattering kernel {^o{£u,£',y) and 
^i{£y^£'^)) via eq. (38) at each point in the computational domain. For 
this reason we tabulate and <hi(e,y,e(,) for each £^ and pair 

on a grid in temperature and ?7e- Because the vector and axial-vector cou¬ 
pling constants for neutrino-electron scattering and the neutrino energy 
grouping differ between Ve-, ^e-, and neutrinos, we construct separate 
tables for each species. Gauss-Legendre quadratures (16-point) are used 
to evaluate the angular integrals over cos 9 for I = 0 and / = 1 in eq. (38). 
At a given temperature/density/composition point the equation of state 
returns rje and we then perform a six-point bivariant interpolation in 
T-r]e space, for the given £u-£l combination, to obtain $o(£!/)£(/) and 
‘hi(e,y, e(,). The integrals over for each energy, which yield and 

j^scatt, then computed using simple trapezoidal rule integration. In 
practice, we use 40 energy groups (nj), 30 temperature points {Np), and 
30 rje points (A^?;). The tables are then I x Uf x n'j x Np x in size, 
with I = 2 (<ho and <I>i), or approximately 50 Megabytes. The source 
and sink at each energy and are then included explicitly 

in a manner analogous to any of the emission or absorption processes. 
Using this method, our calculations including neutrino-electron scatter¬ 
ing, are just 10-15% slower than our calculations ignoring this important 
equilibration process (Thompson, Burrows, & Pinto 2002). 

This calculation of and uses only values of the neutrino 

energy density and flux from the previous timestep and, hence, we intro¬ 
duce an explicit timescale into the energy and electron fraction updates 
returned by our transport algorithm. For this reason, when the scatter¬ 
ing rate is large we may encounter a numerical instability. Because the 
largest scattering rates are encountered when the neutrino phase-space 
distribution function is in local thermodynamical equilibrium and the 
scattering off electrons is unimportant, we simply divide the source and 
sink by a large factor (typically 100 above p = 10^^ g cm“^), thus cir¬ 
cumventing the problem of introducing a short and explicit timescale. 
Again, because fy = at these high densities, this approximation is 
acceptable. 

4.2 Neutrino-Nucleon Scattering 

The kernel for inelastic neutrino-nucleon scattering can be related to 
the non-relativistic structure function employed in the thermalization 
studies of Thompson, Burrows, & Horvath (2000); 


R'^^^{£y,£'^,cos9) = G‘^S{q,(jj)[{l -|- cos0)U^ + (3 — cos0)A^], (62) 
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where S{q,uj) is given in terms of the imaginary part of the polarization 
function, analogous to neutrino-electron scattering (§4.1). 

The neutrino-nucleon scattering kernels, while larger in absolute mag¬ 
nitude than the corresponding neutrino-electron scattering kernels are, 
at most points in energy space and thermodynamic space, much less 
broad. In fact, in thermodynamic regimes most relevant for the forma¬ 
tion of the various species’ spectra, the kernel is quite sharply peaked in 
energy. That is, for a given the distribution of ej^’s is tightly centered 
on £iy because the ratio of the neutrino energy to the nucleon mass is 
small. as a function of cos 6, for neutrino-neutron scattering, is 

shown in Fig. (4) at a representative thermodynamic point. This plot 
is analogous to Fig. (1) for neutrino-electron scattering. By comparing 
the = 19 MeV line with that for = 21 MeV, one sees that the for¬ 
mer is larger and, hence, downscattering is preferred. In addition, the 
overall magnitude is much larger than in the neutrino-electron scattering 
case. Figure (5) shows <ho (eq. 38) for neutrino-neutron scattering as a 
function of for several e,y’s at the same thermodynamic point used in 
Fig. (3), the corresponding figure for neutrino-electron scattering. Note 
that while downscatting is strong for the = 35 MeV kernel, there 
is almost equal upscattering at = 5 MeV. In order to explore the 
effect of this process on the emergent spectra in dynamical simulations 
we must first deal with a technical problem. 

In a typical simulation, we employ 40 energy groups for all neutrino 
species with 1 MeV < < 320 MeV for electron-type neutrinos and 

1 MeV < El/ < 100 MeV for anti-electron and muon neutrinos. The 
grouping is generally logarithmic for the i^eS and linear for i?e and 
Neutrino-electron scattering and neutrino-nucleon scattering are most 
important as thermalization mechanisms at energies below ~ 60 MeV, 
where the phase space distribution of all neutrino species is largest. One 
can see clearly from Fig. (5) that a trapezoidal rule integration of <ho over 

as it appears in eq. (44) and eq. (46) may grossly overestimate 
and la fact, with logarithmic energy grouping one may even 

calculate upscattering when there is none because the energy groups 
become larger with increasing energy. 

For neutrino-electron scattering, we are able to employ a simple trape¬ 
zoidal rule and adequately capture the qualities of the kernel. This im¬ 
plies that we use a linear interpolation for J(, and in each energy bin. 
As Fig. (5) shows, however, in order to get an accurate integral over 
<ho(ei/, e(/), we must do better than simple trapezoidal rule with linear 
interpolation. For neutrino-nucleon scattering, in order to increase the 
accuracy of our scheme without compromising computational efficiency, 
for a given energy grouping we pre-compute a grid of integrals over e(,. 
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We assume that during the dynamical calculation and the computation 
of and both J[, and Hi are proportional to Ael + B over 

an energy interval i+i- Given this assumption and both 

and ^i{eu, el) at a given T and rin,p, we tabulate the following 

integrals; 



where / = 0,1 and each i is an individual energy group, set up at the 
beginning of the calculation. The integrals over el are computed using 
16-point Gauss-Legendre quadrature, each nested with another 16-point 
Gauss-Legendre quadrature for the integral over cos 6 necessary for each 
^i{eu,el). These integrals are tabulated at 30 temperature and 30 rjn^p 
points, analogous to the case for neutrino-electron scattering. For ex¬ 
ample, the integral 



delel^Jl^o{e^,el) 


E 

2 = 1 




i+1 


de',,e',?e[, 


$o(e 




+ 



(63) 


where Ai = {Jl' — Jl'^^')l{el' — ell~^') and Bi = Jl' — Aie]^. The total 
source and sink are given by integrals analogous to eq. (63) with ap¬ 
propriate changes to Ai and Bi, depending on if the term in the source 
or sink is over Jl or Hi- In practice, for a given temperature and den¬ 
sity, we calculate the necessary integrals at the four nearest neighbor 
T — rjn^p points saved in the table and then interpolate the solution using 
a four-point bivariant interpolation scheme. In this scheme, the energy 
integral over the kernel is reproduced extremely well and the primary 
uncertainty in calculating the total source and sink is due to the lin¬ 
ear interpolation of J(, and Hi - the same as in the neutrino-electron 
scattering case. Also interesting is that because the neutrino-nucleon 
scattering kernel is so sharply peaked around e^ and drops off so quickly 
with el, most terms in the sum in eq. (63) are zero. When we fill the 
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table, we note the index i of the lowest and highest intervals that 

contribute significantly (to a part in 10^) to the total integral over 

Typically, only four to five energy groups must be included in the 
final sum. This decreases both the size of the table and the amount 
of time needed to calculate and Although the vector and 

axial-vector couplings for neutrino-neutron scattering are independent 
of neutrino species, we calculate separate tables for t'e, and so as 
to allow for arbitrary energy grouping for each species. 

The Elastic Limit: With the inelastic formalism in hand, it is in¬ 
structive to construct the elastic limit. Note that in the limit of zero 
energy transfer, S{q,uj) —> 5(0) = 2'jTnn,pS{(jj), where Un^p is the neutron 
or proton number density (see eq. 80). The scattering term can then be 
written as 

/* + ! r27T 

- V(l-/0^ d4,'M, 

where M = [{1 + cos6*)E^ + (3 — cos6*)A^]. Using eqs. (36), (41), and 
(42) we find that 

I" d^WJ'SWm - (65) 

where Hin = U^(J(, (iHl) + - nHl/3) and Hout = V‘^{l-Jl- 

-|-3A^ (1 — Jl -|-/3). Integrating over using the delta function, 
we have that 


n scatt 
^(0) 


[/.] 


47rCe2(u2 + 3A2) 


(J. 


fu) + 


I^U2 + 3A2 J 


^n,p^n,p 



fv) 


( 66 ) 


where an,p = {G‘^l'Kc)e^{V‘^ -1-3^2), 8n,p is the scattering asymmetry for 
neutrino-neutron or neutrino-proton scattering, G is the weak coupling 
constant, and c is the speed of light. The result presented in eq. (66) is 
to be compared with the full scattering part of the collision term: 

-t^sfu + ^J ^>(0, ft') Mft') dfl'. (67) 

Performing a first-order Legendre series expansion of the integral elas¬ 
tic scattering source term, and combining the scattering sink, —Ksfu, 


(64) 
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with the scattering source, KgJu + UsSrHv, we obtain eq. (66) and an 
expression for 5n,p- 

Note that taking the zeroth moment of yields zero, as it 

should in the elastic scattering limit. Further, note that the first moment 
is non-zero: 


1 

2 





TTC 


2 



A^) - + 3 ^ 2 ) 


q2 

-+ SA"^) 

TTC 


1 - 


n,p 


^n,p^n,pdiu 


1 ^dn,p 


( 68 ) 


This expression is to be compared with the elastic scattering momentum 
source term on the right-hand side of the first moment of the transport 
equation. In fact, the quantity cTn,pnn,p(l — (5n,p/3) defines the elastic 
transport cross section (see eq. 19). 


5. Dynamic Structure Factors for 
Neutrino—Nucleon Interactions 

An alternate formalism for handling inelastic neutrino-nucleon scat¬ 
tering that more straightforwardly than in §4.2 generalizes to include 
nucleon-nucleon correlations, whether due to fermi blocking or nuclear 
interactions, involves the dynamical structure factor S{q, to). Our discus¬ 
sion here follows closely that found in Burrows and Sawyer (1998,1999). 
Recent explorations into the effects of many-body correlations on neutrino- 
matter opacities at high densities have revealed that for densities above 
10^^ gm cm“^ both the charged-current and the neutral-current inter¬ 
action rates are decreased by a factor of perhaps 2 to 3, depending on 
the density and the equation of state (Burrows & Sawyer 1998,1999; 
Reddy, Prakash, & Lattimer 1998; Yamada 1999). Furthermore, it has 
been shown that the rate of energy transfer due to neutral-current scat¬ 
tering off of nucleons exceeds that due to jz^-electron scattering (Janka 
et al. 1996; Thompson, Burrows, & Horvath 2000). Previously, it had 
been assumed that neutrino-nucleon scattering was elastic (Lamb & 
Pethick 1977). However, these recent reappraisals reveal that the prod¬ 
uct of the underestimated energy transfer per neutrino-nucleon scatter¬ 
ing with cross section exceeds the corresponding quantity for neutrino- 
electron scattering. Since Ve and Pg neutrinos participate in super- 
allowed charged-current absorptions on nucleons, neutrino-nucleon scat¬ 
tering has little effect on their rate of equilibration. However, such 
scattering is important for and equilibration. Since the many- 
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body correlation suppressions appear only above neutrinosphere densi¬ 
ties (~ 10^^ — 10^^ gm cm“^), it is only the kinematic effect, and not the 
interaction effect, that need be considered when studying the emergent 
spectra. In the following we adhere closely to the approach and formal 
development in Burrows and Sawyer (1998,1999). Without interactions, 
the relevant dynamical structure factor, S{q,u}), for neutrino-nucleon 
scattering is simply 


S{q,uj)=2j ■^^J='{\p\){l-J='{\p + q_\))2Tr6{u; + ep-ep+q), (69) 

where .F(|p|) is the nucleon Fermi-Dirac distribution function, Cp is the 
nucleon energy, co is the energy transfer to the medium, and q is the 
momentum transfer. The magnitude of q is related to co and Ei, the 
incident neutrino energy, through the neutrino scattering angle, 6, by 
the expression (see text in §4 after eq. 36), 

q=[El + (El - Lv f - 2Ei{Ei - to) cos 6]^/^ . (70) 


In the elastic limit and ignoring final-state nucleon blocking, S{q,Lo) = 
2Tr6(Lo)nn, the expected result, where is the nucleon’s number density. 

The neutral current scattering rate off of either neutrons or protons 
is (Burrows &: Sawyer 1998), 


d^T 


where 

and 


diod cos 9 

21nc = 


= (47r2)-iG^(Ei - u;f[l - E.(Ei - (u)]Xnc , 

(1 + cos 9)V + (3 — cos 0)^15((7, lo) 


(71) 

(72) 


S{q,co) = 2Imn(°)(l-e"^‘^)“^ (73) 


V and A are the applicable vector and axial-vector coupling terms (see 
§3.4 and §3.5) and f3 = 1/kT. The free polarization function, 11^^^, 
contains the full kinematics of the scattering, as well as blocking due to 
the hnal-state nucleon, and the relevant imaginary part of 11^^) is given 
by: 


Imn*-°^((7,ti;) 


9 

m 

2'Kqf3 


log 


_j_ g—<3++/!/^ 


(74) 



where 


( 75 ) 
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/i is the nucleon chemical potential, and m is the nucleon mass. The 
dynamical structure factor, S{q,u}), contains all of the information nec¬ 
essary to handle angular and energy redistribution due to scattering. 
The corresponding term on the right-hand-side of the transport equa¬ 
tion is: 

S[T^] = {27r)-^Glr J dy.I^c^sF (76) 

where 

Esf = [1 - MEi)]KiEi - - MEi)[l - E^Ei - a;)] (77) 

and p'jy is the final state neutrino momentum. 

In the non-degenerate nucleon limit, eq. (74) can be expanded to 
lowest order in to obtain, using eq. (73), an approximation to the 
dynamical structure factor: 


. n(27rm/3)^/^ _q2 

= — -——e ^+, 

q 

where n is the nucleon number density. This says that for a given 
momentum transfer the dynamical structure factor is approximately a 
Gaussian in uj. 

For charged-current absorption process, e~ -\-p, Imll^^) {q, to) 

is given by a similar expression: 




Imn(°)(g,(j) 


9 

m 

2'n(5q 


log 


_j_ ^ — Q\+Ppn. 


(79) 


Eq. (79) inserted into eq. (73) with a (1 — as is appropriate 

for the charged-current process, substituted for (1 — results in an 

expression that is a bit more general than the one employed to date by 
most practitioners (Bruenn 1985), i.e., S = {Xn — Xp)/{1 — In 

the non-degenerate nucleon limit, the structure factor for the charged- 
current process can be approximated by eq. (78) with n = n^. Note 
that for the structure factor of a charged-current interaction one must 
distinguish between the initial- and the final-state nucleons and, hence, 
between their chemical potentials. To obtain the structure factor for the 
Ue absorption process, one simply permutes pn and pp in eq. (79) and 
substitutes —p for p in the (1 — term. 

In general, the widths of the structure factors are larger than might 
have been expected for scattering off of “heavy” particles. This is be¬ 
cause in the past people thought that the neutrino could lose in 
nucleon scattering an energy equal to only about —El/rrinC^, he. that 
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the fractional energy lost is of order p^/rUnC (~1%). However, this as¬ 
sumes that the nucleons are stationary. In fact, they are thermal and, 
the fractional energy they can in a collision transfer to the neutrino is 
of order pn/mnC. Since the nucleons have such a large mass, if they 
and the neutrino have the same energy, pn/rrinC is much larger than 
Vvl'f^nC, at incident neutrino energies of 10-30 MeV by as much as an 
order of magnitude. The formalism above incorporates the kinematics 
of such a collision, a realistic Fermi-Dirac energy distribution for the 
nucleons, and final-state nucleon blocking. The upshot is broad distri¬ 
butions. Including many-body effects further flattens and broadens the 
distribution, while lowering the central values of da/duj, as well as the 
total integral over to (Burrows & Sawyer 1998,1999). 

5.1 Aside: Static Structure Factors 

In the limit of heavy nucleons, when we perform the integration in eq. 
(69) over a range of los and evaluated the inner factors at w = 0, we can 
express this limit as, 

(27r)“^5((7,a;) —> (27r)“^(5(w) J dev'S{q,uj') = 6{co)S{q), (80) 

where S{q) is the static structure factor. At the high densities and 
temperatures achieved in the supernova context, the co = 0 (elastic) limit 
is not particularly accurate (Burrows & Sawyer 1998; Reddy, Prakash, 
& Lattimer 1998). 

S{q) is merely the Fourier transform of the thermally-averaged density- 
density correlation function. This is the classic result that scattering off 
of a medium is in reality scattering off of the fluctuations in that medium. 
Also of interest is the long-wavelength limit, (7 —> 0, justified when the 
neutrino wavelength is much bigger than the interparticle separation. 
Statistical mechanics provides two useful and equivalent expressions for 
the long-wavelength limit, S'(O), the first (Landau & Lifschitz 1969) in 
terms of the isothermal compressibility of the medium Kt (= _ |^)^ 

5(0) = = n^, (81) 

Ao 

where Kq is the ideal gas compressibility and h is the average nucleon 
density, and the second in terms of the derivative of the density with 
respect to the chemical potential of the nucleons, p, 




( 82 ) 
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In the ideal gas limit of no correlation between particles, eqs. (81) and 
(82) show that S{q) is simply equal to the number density, h, as ex¬ 
pected from eq. (69), without blocking. Eq. (81) reveals that if Kp 
is small because the matter is stiff, in the long-wavelength limit the 
neutrino-matter cross sections are suppressed. In general, when we ac¬ 
knowledge that the neutrino-matter interaction has axial-vector current 
and nucleon isospin terms, we require separate correlation functions for 
the neutron and the proton, as well as for spin correlations. These de¬ 
pend upon susceptibilities that are different from the compressibility, 
but we find suppression in these terms as well (Burrows & Sawyer 1998; 
Reddy, Prakash, & Lattimer 1998). 

Eq. (82), equivalent to eq. (81) by a thermodynamic identity, is a pow¬ 
erful result of great generality. In standard approximation schemes for 
the many-body problem, the distribution function for a nucleon species 
is given by the Fermi-Dirac distribution in which the chemical potential 
fi is replaced by — where v{n) is the average energy of interaction 
of the nucleon with the other nucleons and is a function of the density. 
Thus, the density is given implicitly by 



The expression (83) holds in the Hartree approximation; it holds in 
approaches that introduce mean meson fields instead of nuclear poten¬ 
tials; it holds in the Landau Fermi liquid theory (FLT), subject to the 
proviso that we use only results in which the derivative of the potential 
V (with respect to the n) enters; and it holds in approaches using the 
Skyrme potential. 

Differentiating (83), we can solve for and S'(O), 

^( 0 ) = = h{f,)[i + HR)^r\ ( 84 ) 

where 


Hr) = 2 y 


(Pp gd\p'^/{'2ni)-n+v] 

(27r)3 [1 -I- e/^[p^/(2"i)-/^+d]2 


= 2 


d^p 


J^{p){l-J^ip)) 


(85) 

and J-{p) is the Fermi-Dirac function, but with the chemical potential 
displaced by u. If we regard particle densities as inputs to our calcu¬ 
lations, then the displacement of the chemical potential by the nuclear 
potential is irrelevant, since the same difference, p- — v, enters the cal¬ 
culation of the density in terms of the chemical potential. Thus, the 
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numerator of (84) contains no more than the familiar Pauli blocking ef¬ 
fects (for the case q = 0); the denominator contains all of the effect of 
the interactions. 

As an example, consider a two-nucleon potential V(r). In the Hartree 
approximation, the average potential seen by a single nucleon is given 
by u = nil, where U = f d^xV{x), and (84) becomes 

S{0)=h{iJ.)[l + h{fj,)U]~\ (86) 

the potential providing an enhancement, if negative, and a suppression, 
if positive. The latter is the case for high-density nuclear matter. 

5.2 Procedure for Calculating i/—nucleon 

Structure Functions for Neutral-Current 
Scattering Including Interactions 

Now, taking only the neutron part of the vector-current coupling, the 
differential neutrino-nucleon scattering rate is given by. 


d^r 

duj d cos 9 


- U{E2)] {l + COSe){C^YSnn{q,00 


-F(3 - cose)g\[S^p{q,uj) + - 25^„(g,w)], (87) 


pn V 


where E 2 =Ei — ui. 

The structure functions, S (Fermi) and (Gamow-Teller; axial), are 
elements of separate 2x2 symmetric matrices. For the vector dynamic 
structure function, S, we have 


S{q,uj) 


f Spp{q,io) Spn{q,uj)\ 
Snn{,Qi^') ) 


The structure function matrix is given by. 


S'(g, w) = 2Im n® (g, w)[l — u(( 7 )n® (( 7 ,ti;)] ^ (1 — e 




( 88 ) 


where 


nW(g,u;) 


/n®(g,u;) 0 \ 

Iv 0 n®(g,u;)y 


and IIp^^ and Ili^^ are given by the polarization function and evaluated 
with the proton and neutron chemical potentials, respectively (eq. 74). 
The potential matrix is. 
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^ ^ / vi+ V2 + 4:TTe^{q‘^q^p) ^ -ui - ^2 \ 

V Vi — V2 f 1 + f2 / 

where the v's were defined in terms of Fermi liquid parameters and the 
term containing q^p is the Thomas-Fermi screened Coulomb potential 
{Qtf — 4e^7r^/^(3hp)^/^). Following Burrows and Sawyer (1998,1999) 
and for simplicity, we use here FLT and Landau parameters, in lieu of a 
more developed nuclear interaction model. 

In a real calculation, in all the kinematic expressions the nucleon mass 
(m) is to be replaced by m*. Unfortunately, the relation of Landau 
parameters to experimental results depends upon the effective mass in 
model-dependent ways. Taking m* = 0.75mn as our fiducial value for 
the effective mass, we use parameters from Backman, Brown, &: Niskanen 
(1985) and Brown and Rho (1981): Fq = —0.28; Fq = 0.95; Go = 0; Gq = 
1.7 and A = 2.63 x 10“®MeV“^, obtaining, 

vi = -7.4 X 10-® MeV“2 
V2 = 2.5 X 10"® MeV^ 

Vs = 0 

Vi = 4.5 X 10-5 MeV-2. (89) 

For other values of the effective mass, we keep these potentials at the 
same value, which is to say we assume that the Landau parameters are 
proportional to m*/m. 

The form for the Gamow-Teller matrix, S^{q,u), is the same as that 
for S, except that the potential matrix is replaced by v^ 

yA ^ ( Vs + Vi U3 - U4 A 
\Vs-Vi Vs + Vi ) ' 

Taking the matrix inverses leads to the following forms for the com¬ 
binations of structure functions that appear in (87) 

Snn{q,io) = 2Im[nWzi-i](l - (90) 

where 

Dy = l-{vi+ U 2 )nW - (Ul - U2)2nWn(0)Qv'-l . (91) 

Qv is given by the expression: 

Qv = 1- 47re^(g^ qpp)~^U^^^ - {vi ■U2)n® . 


(92) 
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If, as in (89), we take V 3 = 0, we obtain the simple result for the 
axial-current terms, 


SAiq,u}) = 2 Im 


{q, uj) + n))''’ (g, uj) 


( 0 )/ 


1 - U4 (g, Oj) + (g, uj)] 


( 0 )/ 


(1 - . (93) 


For the Fermi term, since Cy'^ = 1/2 — 2sin^0ru ~ 0; we drop the 
proton structure function in (87). Furthermore, we use the potential 
parameters given in eq. (89), and in eq. (91) we drop the third term. 
This term would have been significant had it not been for the Coulomb 
term in the denominator, an illustration of the importance of the explicit 
inclusion of Coulomb forces, even for the neutron density correlations. 
Since the ViS are all real, we obtain for the structure factors used in (87), 



SF{q,uj) = 2ImnW(i _ 

(94) 

where 




Cy = (1 - UFRenW )2 + 

(95) 

and 




SA{q,u;) = 2[lmnf((7,(u) + ImnW(g,u;)] (1 - 

where 

Ca = ^Ai + 1^.42 • 

and Ca 2 are given by the expressions; 

= [1 - UGT(Renf(g,(u) + RenW(g,cu)) 


and 


Ca 2 = Vqt [imR® {q, cv) + Imn^°) {q, uj] 


(96) 

(97) 


(98) 

(99) 


The F in SF{q, uj) and the A in SA{q, oj) stand for Fermi and Gamow- 
Teller (axial) and vp and vgt equal (ui + V 2 ) and U 4 , respectively, in 
Fermi Liquid Theory. SA{q,^) in eq. (96 is now the entire axial term 
in eq. (87). is the correction factor due to many-body effects for a 
given momentum transfer (or scattering angle) and energy transfer. A 
similar procedure is employed for calculating the many-body corrections 
to the charged-current rates (Burrows & Sawyer 1999). 
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5.3 Collective Excitations of the Medium 

Following Burrows and Sawyer (1998,1999), we note that for most 
regions of phase space, Cy and in eqs. (95) and (97) are greater than 
one and represent suppression in the scattering rates. Their effects on 
the integrals over uj and 0 are always suppressive. However, the terms 
containing the real parts have roots; these roots represent collective exci¬ 
tations. For the Fermi term, zero sound in the medium can be generated 
if the scattering has a {co, q) pair that satisfies the mode’s dispersion rela¬ 
tion, i.e., if it hits the resonance. Similarly, for the Gamow-Teller term, 
spin waves in the protons and the neutrons (related by a set phase) can 
be generated. These modes are the traveling-mode equivalents of the 
Gamow-Teller resonance in nuclei (a standing wave). The zero sound 
of the Fermi part is analogous to the Giant-Dipole resonance in nu¬ 
clei. The resonances increase the structure function when the scattering 
transfer ratio, uj/q, equals the ratio of the collective excitation’s angu¬ 
lar frequency and wavenumber. For a given scattering angle, one can 
plot the differential cross section in lv and cos 9 as a function of eo/q 
to see the resonances. In Figure (6), we display this for five different 
angles between 15° and 180°, an incident neutrino energy of 20 MeV, a 
temperature of 5 MeV, a density of 3 x 10^^ g cm“^, and an electron 
fraction, T)., of 0.3. We see in Figure 6 that the resonances in both the 
forward and the backward directions line up at the same values of oj/q, 
as expected for a collective mode, and we can straightforwardly calcu¬ 
late the mode’s dispersion relation. This is akin to the Cerenkov effect. 
Note that the Gamow-Teller term dominates the Fermi term, so that in 
Figure 6 we are really seeing the spin waves related to the Gamow-Teller 
resonance. However, the dispersion relations for zero sound and these 
spin waves are generally similar. In fact, recalling the classic result of 
Fetter and Walecka (1971) that in the weak-coupling limit, the speed of 
zero sound in a degenerate system is ~ Vfermi, where Vfermi is the Fermi 
velocity, and recalling that for nucleons in nuclei Vfermi is ~ 0.3c, the 
calculated resonance value of oj/q is not unexpected. In Figure 7, we 
plot the Gamow-Teller structure function versus oj/q for various values 
of a;, m*, and two values of the density. At m* = for each value 
of the density we obtain a sharp resonance, but at two different speeds, 
reflecting the crude p^/^-dependence expected for Vfermi- For a given 
density, the mode speed is seen in Figure 7 to be inversely proportional 
to the effective mass. The width of the resonance is determined by the 
magnitude of the imaginary part of the polarization function. 
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6 . 


e+e Annihilation 


Ignoring phase space blocking of neutrinos in the final state and taking 
the relativistic limit {me —> 0), the total electron-positron annihilation 
rate into neutrino-antineutrino pairs can be written in terms of the 
electron and positron phase space densities (Dicus 1972); 

= Ki f [ J='e-J='e+{4-4+ +4-4+) dee- dee+ , 


me+ 


Tic 


( 100 ) 

where Ki = (l/187r^)ccro(C'y + C*^)- Again, Cy = 1/2 + 2sm^0w for 
electron types, Cy = —1/2 + 2sin^ 9w for mu and tau types, and = 

(1/2)^. Rewriting eq. (100) in terms of the Fermi integral Fn{ri), we 
obtain: 

Q,^,^=K,{kT)^^j [FMFs{-Ve) + F,i-r^e)Fs{rie)] , 

( 101 ) 

where rje = f+e/kT and 


roo 

Fniv) = / - dx . 

^ " Jo + 1 


Integrating eq. (100), we obtain 


~ 9.7615 X 10 


24 


kT 


MeV 


-3^-1 


f{r]e) ergs cm '’s 


where 


fide) = 


Fiirie)Fo{-rie) + Fi{-rie)Fo{'qe) 
2F4(0)F3(0) 


For Vui+a and Vt-Vt production combined, 


Q. 


~ 4.1724 X 10^^ 


kT 

MeV 


fi+ie) ergs cm '^s 


- 3,-1 


( 102 ) 


(103) 


(104) 


(105) 


One can easily derive the spectrum of the total radiated neutrino 
energy (et) by inserting a delta function (/ 5{eT — Se- — Se+)deT) into 
eq. (100). Recall that the total energy of the neutrinos in the hnal state 
is equal to the sum of the electron and positron energies in the initial 
state. Integrating first over £^+ to annihilate the delta function and then 
over Eg- to leave a function of et, one obtains; 



^e— ] dEe— . 
(106) 
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The numerical evalution of eq. (106) is straightforward. The average of 
et is equal to; 


(er) = ( 


FiiVe) 

FsiVe) 


+ 


F4{-r]e) \ 

Fsi-Ve)^ 


(107) 


which near r/e ~ 0 is ~ 8 T and for r/e >> 1 is ~ 4T(1 + rje/B). 

However, while the total energy loss rate (eq. 103) and the spectrum 
of Et pose no great mathematical problems, the production spectrum of 
an individual neutrino is not so easily reduced to a simple integral or to 
an analytic expression. This is due primarily to the awkward integration 
of the angular phase space terms, while subject to the momentum con¬ 
servation delta function, and to the explicit dependence of the matrix 
elements on the electron/neutrino angles. /,From Dicus (1972), averag¬ 
ing over initial states and summing over final states, the matrix element 
for the e'’'e“ ^ i'i> process in the me = 0 limit is; 


^ ^ |>1 P = leG^ [{Cv + Ca^P • qp p' • Hi. + {Cv - Ca^P • p' • qp] , 

(108) 

where p and p' are the four-momenta of the electron and positron, re¬ 
spectively, and Qi, and qp are the four-momenta of the neutrino and 
antineutrino, respectively. Using the formalism of Bruenn (1985) and 
Fermi’s Golden rule, expanding the production kernel in the traditional 
truncated Legendre series, performing the trivial angular integrals, tak¬ 
ing the non-trivial angular integrals from Bruenn (1985), and ignoring 
final-state neutrino blocking, we obtain for the single-neutrino source 
spectrum due to e+e“ annihilation; 


where 

and 


dQ 

dE„ 


87r^ 3 

(27r;ic)6 



dEpEl^l{Eu,Ep) , 


(109) 


TT 7o 


dEe-F 5- {Ee-]F =+ [Eu -\- Ep — Eg-] Hq{e,^, EpEe-) , 

( 110 ) 


Ho{Ei^, Ep, eg- ) — {Cv+c a)‘^ Jq{£vCPi ^e-)F{Cv — Ca)^ Jq^ {e^, Ep, Cg- ) . 

(Ill) 

The Jqs in eq. (Ill) come from the more obdurate angular integrals 
required by the dot products in eq. (108) and the momentum delta 
function and have the symmetry; 


do {Ei/, El/, Ee—') Jq (Ei/, El/, Ee—} . 


( 112 ) 
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/.From eqs. (109) and (111), we see that the differences between the 
spectra of the Vf. and neutrinos flow solely from their correspondingly 
different values of {Cy + Ca)'^ and {Cy — Ca)'^- One can use 4-point 
Gauss-Legendre integration to calculate eq. (110) and 16-point Gauss- 
Laguerre integration to calculate eq. (109). 

At small r]e, the e'’'e“ annihilation spectra and total energy loss rates 
for the i^e and Pe neutrinos are similar, as are the average emitted i^e 
and Fe neutrino energies. However, as ije increases, both the total en¬ 
ergy radiated in Fg neutrinos and the average Pg energy start to lag the 
corresponding quantities for the i^e neutrinos. This is true despite the 
fact that the total number of i^e and Fg neutrinos radiated is the same. 
If final-state blocking is ignored, {£i)/T is a function of r]e alone, be¬ 
coming linear with rje at high r]e and one half of eq. (107) (~4.0) at 
low r]e- Note also that {e,^^)jT and are closer to one another 

than are {s^^)/T and {£p^)fT. The individual production spectra vary 
in peak strength, in peak energy, and in low-energy shape, but they are 
quite similar on the high-energy tail. Due to the parity-violating matrix 
element for the e+e“ ^ vu process and the fact that r]e is positive, the 
antineutrino spectra of all species are softer than the neutrino spectra. 
The pair sums of the integrals under these curves are given by eqs. (103) 
and (105). For r]e = 0, 50% of the pair energy emission of electron types 
is in i?e neutrinos, but at r]e = 10 only 42% of this total energy is in Fg 
neutrinos. However, at rje = 10, the neutrinos still constitute 48.5% 
of the pair emission. These differences reflect differences in the 

corresponding coupling constants Cy and Ca- 


7. UiUi Annihilation 

In the limit of high temperatures and ignoring electron phase space 
blocking, the yii>i annihilation rate into e'’'e“ pairs can be written (Janka 
1991): 



where Jiy is the zeroth moment of the radiation field, £i, is the neutrino 
energy, Ki is defined as before {i.e., Ki = (l/lSvr^)ccro(C'y -|- G^)), and 


^ +PuiPPi + -{1 - Pui)il - PPi) 

(114) 

where the flux factor and the Eddington factor = 

{Pu ) — Pvl^v- Eq. (113) can be rewritten in terms of the invariant 
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distribution functions 



el, + etet) dep, . 


(115) 


Note that when the radiation field is isotropic = 1) and when 
rje = 0 the total rate for e'’'e“ annihilation given in eq. (100) equals 
that for r'ji/j annihilation given in eq. (115), as expected. Buras et al. 
(2002) have addressed the related and interesting process of —> i^jUj- 

We refer to that paper for a discussion of the relevance and rates of this 
process. 

8. Nucleon—Nucleon Bremsstrahlung 

A production process for neutrino/anti-neutrino pairs that has re¬ 
cently received attention in the supernova context is neutral-current 
nucleon-nucleon bremsstrahlung (ni-|-n 2 ^ U 3 -|-n 4 -|-z^z/). It importance 
in the cooling of old neutron stars, for which the nucleons are quite de¬ 
generate, has been recognized for years (Flowers 1975), but only in the 
last few years has it been studied for its potential importance in the 
quasi-degenerate to non-degenerate atmospheres of protoneutron stars 
and supernovae (Suzuki 1993; Hannestad & Raffelt 1998; Burrows et al. 
2000; Thompson, Burrows, & Horvath 2000). Neutron-neutron, proton- 
proton, and neutron-proton bremsstrahlung are all important, with the 
latter the most important for symmetric matter. As a source of i^e and 
i?e neutrinos, nucleon-nucleon bremsstrahlung can not compete with the 
charged-current capture processes. However, for a range of temperatures 
and densities realized in supernova cores, it may compete with e^e“ an¬ 
nihilation as a source for i/t, and i?r neutrinos (“t'^”s). The 

major obstacles to obtaining accurate estimates of the emissivity of this 
process are our poor knowledge of the nucleon-nucleon potential, of the 
degree of suitability of the Born Approximation, and of the magnitude 
of many-body effects (Hannestad & Rafflet 1998; Raffelt &: Seckel 1998; 
Brinkmann & Turner 1988). Since the nucleons in protoneutron star at¬ 
mospheres are not degenerate, we present here a calculation of the total 
and differential emissivities of this process in that limit and assume a 
one-pion exchange (OPE) potential model to calculate the nuclear ma¬ 
trix element. For the corresponding calculation for arbitrary nucleon 
degeneracy, the reader is referred to Thompson, Burrows, & Horvath 
(2000). The formalism we employ has been heavily influenced by those 
of Brinkman and Turner (1988) and Hannestad and Raffelt (1998), to 
which the reader is referred for details and further explanations. 
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Our focus is on obtaining a useful single-neutrino final-state emission 
(source) spectrum, as well as a final-state pair energy spectrum and the 
total emission rate. For this, we start with Fermi’s Golden Rule for the 
total rate per neutrino species: 



d^Pi 1 d^qp 

(27r)3- (27r)^2u!,y (27r)^2wp 




where 


^brems = (116) 

(5^(P) is four-momentum conservation delta function, uj is the energy 
of the final-state neutrino pair, {ujy,qu) and {ujp,qp) are the energy and 
momentum of the neutrino and anti-neutrino, respectively, and p) is 
the momentum of nucleon i. Final-state neutrino and anti-neutrino 
blocking have been dropped. 

The necessary ingredients for the integration of eq. (116) are the 
matrix element for the interaction and a workable procedure for handling 
the phase space terms, constrained by the conservation laws. We follow 
Brinkman and Turner (1988) for both of these elements. In particular, 
we assume for the n + n^n + n + uP process that the matrix element 
is: 




64 


G if/mn) Pa 




+ ml 


f + 


= A 


OJyOJp 
2 ’ 


(117) 


where the 4 in the denominator accounts for the spin average for identical 
nucleons, G is the weak coupling constant, / (~ 1.0) is the pion-nucleon 
coupling constant, qa is the axial-vector coupling constant, the term in 
brackets is from the OPE propagator plus exchange and cross terms, k is 
the nucleon momentum transfer, and is the pion mass. In eq. (117), 
we have dropped iji, ■ k terms from the weak part of the total matrix 
element. To further simplify the calculation, we set the “propagator” 
term equal to a constant a number of order unity, and absorb into C 
all interaction ambiguities. 

Recently, Hanhart, Phillips, & Reddy (2001) have addressed these mo¬ 
mentum terms in the context of axion emission and production in 
supernovae. In an effort to make contact with the approximation to the 
matrix element we present here, they plot C as a function of average rel¬ 
ative thermal nucleon momentum (p; Phillips, private communication). 
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The function peaks for C,{p) between 150 — 200 MeV at C — 0.47. At 
p = 50 MeV C ~ 0.08 and at p = 500 MeV ( ~ 0.27. We are most inter¬ 
ested in the region around the neutrinospheres, where the emergent 
spectrum might be most affected by nucleon-nucleon bremsstrahlung. 
Mass densities and temperatures in this region might be 10^^ — 10^^ g 
cm“^ and 5 — 10 MeV, respectively. We estimate p in this regime to 
be ~ 175 MeV and take C = 0-5 for all thermodynamical points. The 
constant A in eq. (117) remains. 

Inserting a / 5{uj — lOy — ujp)d(jj by the neutrino phase space terms 
times LoiOyLOyliJ^ and integrating over wp yields: 

r LOyUJy cPqp 1 /■“ Ulljui -LVuf ^ 

J ^ uj'^ { 27 r)^ 2 (jj^ { 27 r)^ 2 (jjp ( 27 r)^ Jo Jo ^, 

(118) 

where again to equals {uju top). If we integrate over we can derive 
the CO spectrum. A further integration over to will result in the total 
volumetric energy emission rate. If we delay such an integration, after 
the nucleon phase space sector has been reduced to a function of co and 
if we multiply eq. (116) and/or eq. (118) by coy/co, an integration over co 
from cOu to infinity will leave the emission spectrum for the single final- 
state neutrino. This is of central use in multi-energy group transport 
calculations and with this differential emissivity and Kirchhoff’s Law 
(§ 2 ) we can derive an absorptive opacity. 

Whatever our final goal, we need to reduce the nucleon phase space 
integrals and to do this we use the coordinates and approach of Brinkman 
and Turner (1988). We define new momenta; = (pi -|-p 2 )/ 2 , P- = 
(pi — P2)/2, P 3 c = Po — P+-, and p^c = Pa — P+-, where nucleons 1 and 2 
are in the initial state. Useful direction cosines are 71 = ■p-/\p+\\p-\ 

and 7c = p+ • Pi,c/\p+\\pzc\- Defining Ui = pf/2mT and using energy and 
momentum conservation, we can show that: 


d^Pid^P2 

= 8d^p+d^p- 


CO 

= 2T{u--usc) 


Ul^2 

= n+-|-tt_ ± 2 (n+n_)^/^ 7 i 


'W3,4 

= n+-Ftt 3 c± 2 (u+n 3 c)^/^ 7 c. 

(119) 


In the non-degenerate limit, the — ^a) term reduces 

to where y is the nucleon degeneracy factor. Using eq. 

(119), we see that the quantity (tt+ U-) is independent of both 71 
and 7 c. This is a great simplification and makes the angle integrations 
trivial. Annihilating d^p^ with the momentum delta function in eq. 
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(116), noting that pfdp = ^ dui^ pairing the remaining energy 

delta function with tt_, and integrating from 0 to oo, we obtain: 

dQnb = ^ — 5 / e~'"{x^+xuj/T)^/‘^dx (hj. 
2* X 3 X LJo i 

( 120 ) 

The variable x over which we are integrating in eq. (120) is equal to 
2u^c- That integral is analytic and yields; 



e ^{x^ + xui/T^/'^dx = 


r]e'^Ki{r]) , 


( 121 ) 


where Ki is the standard modified Bessel function of imaginary argu¬ 
ment, related to the Hankel functions, and rj = ujjTT. Hence, the lo 
spectrum is given by: 


dQnb 

duj 




( 122 ) 


It can easily be shown that (cu) = 4.364T. Integrating eq. (120) over 
u) and using the thermodynamic identity in the non-degenerate limit: 



where is the density of neutrons (in this case), we derive for the total 
neutron-neutron bremsstrahlung emissivity of a single neutrino pair; 


Qnb = 1-04 X 10^°C(-^n/5i4)^(7i-^)^'^ergscm ^s ^ 

Ivlev 


where pn is the mass density in units of 10^^ gm cm“^ and Xn is the 
neutron mass fraction. Interestingly, this is within 30% of the result in 
Suzuki (1993), even though he has substituted, without much justifica¬ 
tion, (l-|-ti;/2T) for the integral in eq. (120). ([1-|-(7r77/2)^/^] is a better 
approximation.) The proton-proton and neutron-proton processes can 
be handled similarly and the total bremsstrahlung rate is then obtained 
by substituting -|- X^ -\- ‘^X^Xp for X"^ in eq. (124) (Brinkmann and 
Turner 1988). At Xn = 0.7, Xp = 0.3, p = 10^^ gm cm“^, and T = 
10 MeV, and taking the ratio of augmented eq. (124) to eq. (105), we 
obtain the promising ratio of ~ 5((. Setting the correction factor C, equal 
to ~ 0.5 (Hanhart, Phillips, and Reddy 2001), we find that near and 
just deeper than the Up neutrinosphere, bremsstrahlung is larger than 
classical pair production. 
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If in eq. (118) we do not integrate over uj^, but at the end of the 
calculation we integrate over lo from lOu to oo, after some manipulation 
we obtain the single neutrino emissivity spectrum; 

= ( 125 ) 

= ( 126 ) 

where r]i, = (jj^/2T, C is the normalization constant equal to 
(= 0.564), and for the second expression we have used the integral rep¬ 
resentation of Ki{rf) and reversed the order of integration. In eq. (126), 
Qnb is the emissivity for the pair. 

Eq. (126) is the approximate neutrino emission spectrum due to 

nucleon-nucleon bremsstrahlung. A useful fit to eq. (126), good to 

better than 3% over the full range of important values of r]^, is: 

dQ'nb ^ ( 127 ) 

diOu T V T / 

Thompson, Burrows, and Horvath (2000) should be consulted for a de¬ 
tailed discussion of nucleon-nucleon bremsstrahlung for arbitrary nu¬ 
cleon degeneracy. 

9. Conclusion 

The processes that have been described above are essential elements 
of the neutrino-driven supernova explosion mechanism. Coupling these 
with radiation-hydrodynamics codes, an equation of state, beta-decay 
and electron capture microphysics, and nuclear rates, one explores the 
viability of various scenarios for the explosion of the cores of massive 
stars (Liebendorfer et al. 2001ab; Rampp & Janka 2000). Recently, 
Thompson, Burrows, and Pinto (2002) have incorporated this neutrino 
microphysics into simulations of ID (spherical) core collapse and have in¬ 
vestigated the effects on the dynamics, luminosities, and emergent spec¬ 
tra of weak magnetism/recoil, nucleon-nucleon bremsstrahlung, inelastic 
neutrino-electron scattering, and a host of the cross section corrections 
described above. Figures 8 and 9 depict some of the resultant luminos¬ 
ity spectra and their temporal evolution for a representative simulation. 
The character of the spectra reflect the opacities and sources described in 
this paper. In particular, the energy hardness hierarchy from (softer) 
to (harder) neutrinos is clearly demonstrated on these plots, as is 
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the distinction between the Vf. pre-breakout and post-breakout spectra 
(Fig. 8). (See the figure captions for further details.) 

To date, none of the detailed ID simulations that have been performed 
explodes and it may be that multi-dimensional effects play a pivotal role 
in the explosion mechanism (Herant et ah 1994; Burrows, Hayes, &: 
Fryxell 1995; Fryer et ah 1999; Janka & Muller 1996). Be that as it 
may, an understanding of neutrino-matter interactions remains central 
to unraveling one of the key mysteries of the nuclear universe in which 
we live. 
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Figure 2. For the same thermodynamic point as used for Fig. (1), the scattering 
kernel thick solid lines) for electron scattering as a function of cos 6, for £i, = 

20 MeV and — 2, 10, and 16 MeV. Short dashed lines show the first-order Legendre 
series expansion approximation to which is linear in cosO; 7?°“* ~ (l/2)<Fo + 

(3/2)'Fi cos9. The long dashed line shows the improvement in going to second order 
in COS0 by taking ~ (l/2)<l?o + (3/2)'Fi cos6 (5/2)<F2(l/2)(3cos^ 6 — 1). 
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Figure 3. The I — 0 term in the Legendre expansion of the electron scattering 
kernel, il?o(ei^,eL) (eq. 38), for T = 6 MeV and r/e = 20 as a function of e'^ for = 5, 
10, 15, 20, 25, and 35 MeV. Note that for any e^, the neutrino is predominantly 
downscattered. The magnitude of <E>o(e,/,eG and sign of (w) are to be compared with 
those in Fig. (5). 
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Figure 4- The scattering kernel R°'^^£v,cos6) for i^e—neutron scattering as a 
function of cos 6 for = 20 MeV and = 18, 19, 21, and 22 MeV, at a representative 
thermodynamic point (T = 5 MeV, p = 10^® g cm“®, = 0.5). Note that although 

the absolute value of the energy transfer (|e^ — e),j) is the same for both e'^, = 19 
MeV and = 21, the absolute value of i?°“*(20,19, cos 0) is greater than that of 
i?°“*(20, 21, COS0), reflecting the fact that at this temperature the incoming neutrino 
is more likely to downscatter than upscatter. 
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Figure 5. The 1 = 0 term in the Legendre expansion of the neutrino-nucleon scat¬ 
tering kernel, <l>o(e,/,eO (®<1- 38), for T = 6 MeV and r;„ = —2 as a function of 
for Ev = 5, 10, 15, 20, 25, and 35 MeV. Note that for = 5 MeV the neutrino is 
predominantly upscattered, while for Ei, = 35 MeV the neutrino is predominantly 
downscattered. The magnitude of <l>o(ei/, eL) a-nd sign of (ui) are to be compared with 
those in Fig. (3). 
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u/q 


Figure 6. Logio of the doubly-differential cross section for neutral-current 
neutrino-nucleon scattering versus uj/q for scattering angles 15°, 45°, 90°, 135°, and 
180° . The calculations were performed at a temperature of 5 MeV, a Ye of 0.3, a p of 
3 X 10^“^ g cm“®, and an incident neutrino energy of 20 MeV. The default potentials 
{vgt = 4;.5 X 10“® and vf = 1.76 x 10“®) and effective mass (m* = 0.75 m„) were 
employed. The differential cross section is divided by the total scattering cross section 
(cti) in the non-interacting, no-nucleon-blocking, a; = 0 limit. (Figure taken from 
Burrows and Sawyer 1998.) 
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Figure 7. Logio of the Gamow-Teller structure function versus uj/q for an incident 
neutrino energy of 20 MeV, energy transfers, u), of 6, 8, and 10 MeV, two values of 
the effective mass (m* = [0.75m„, l.Omn]) and two values of the density (p = 3 x 10^^ 
and 10^® g cm“®). A temperature of 5 MeV and a Ye of 0.3 were used, as was the 
default vgt (= 4.5 x 10“®). (Figure taken from Burrows and Sawyer 1998.) 
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Figure 8. A collection of emergent spectra at various times during the core- 
collapse, bounce, and shock-stagnation phases of the core of an 11 Mq progenitor. 
The luminosity spectrum (logarithm base ten) is in units of 10®"* ergs s“^ MeV”*^ 
and the neutrino energy (abscissa) is in units of MeV. The dashed curves cover the 
collapse phase (of duration ^200 milliseconds) until just before the peak luminosity 
around shock breakout is achieved and the solid curves are for the subsequent cooling 
and deleptonization phases after the peak. The last curve is at ~110 milliseconds 
after bounce. The lowest dashed curve is at a time early during collapse. Note the 
relative softness of the spectrum then. As the figure shows, the transition from the 
dashed to the solid curve happens close to the time when the spectrum is hardest. 
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Figure 9. This figure shows the Ve (thin lines) and (thick lines) emergent 

luminosity spectra for the 11 Mq progenitor evolution depicted in Fig. 8. The 
luminosity spectra (logarithm base ten) are in units of 10®"^ ergs s“^ MeV“^ and 
the neutrino energy (abscissa) is in units of MeV. There is no appreciable flux prior 
to shock breakout for these species. To avoid clutter, we here depict only a few 
spectra to ~50 milliseconds after bounce. (These curves represent the sum of the 
Vfj,, v-r, and Vt luminosity spectra.) However, the t'e spectra are shown until about 
110 milliseconds after bounce. During the phases shown, both sets of luminosities are 
always increasing. Note that the Vfj, spectra are significantly harder than either the 
Ve or the Ve spectra. This is a consequence of of the fact that the i^^s do not have 
appreciable charged-current cross sections (eqs. 10 and 11), enabling one to probe 
more deeply into the hot core with these species. 




